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ABSTRACT 

Scattering  of  body  waves,  longitudinal  and  shear,  by  a  normal  surface  breaking  crack  has 
been  studied  in  this  paper.  A  hybrid  method  that  combines  the  integral  representation  of 
the  scattered  field  with  a  finite  element  discretization  of  the  near-field  has  been  used  to 
study  the  crack  opening  displacement  (COD)  and  the  free-surface  displacement  (SD).  It  is 
shown  that  both  COD  and  SD  are  insensitive  to  the  crack-tip  field.  In  addition,  at  low 
frequencies  the  COD  is  found  to  be  uniform  over  most  of  the  crack  length  and  thus  can 
be  approximated  by  a  constant,  which  may  be  estimated  from  far-field  observations^^ 
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INTRODUCTION 

Problems  of  elastic  wave  scattering  by  surface-breaking  and  near-surface  cracks  are  of 
considerable  current  interest  for  ultrasonic  nondestructive  evaluation.  Ultrasonic  scattering 
by  planar  cracks  near  or  at  the  free  surface  of  a  semi-infinite  elastic  isotropic 

homogeneous  medium  has  been  studied  by  many  authors  (Shah  et  al..  1985.  1986; 

Achenbach  ei  at.,  1984;  Van  der  Hijden  and  Neerhoff.  1984;  Shah  et  at.,  1987;  Zhang 

and  Achenbach.  1988a.b). 

In  this  paper  we  analyze  the  problem  of  body  wave  (plane  longitudinal  and  shearl 

scattering  tjy  a  surface-breaking  crack  with  particular  emphasis  on  the  crack  opening 
displacement  (COD)  and  the  free  surface  displacement  (SD)  near  the  mouth  of  the  crack. 

Our  object  here  is  to  assess  the  dependence  of  these  quantities  on  the  representation  of 

the  near-the-crack-tip  displacement  field  in  a  finite  element  formulation.  It  is  now  well 
established  that  in  order  to  capture  the  crack-tip  stress  singularily  in  a  finite  element 

formulation  it  is  necessary  to  use  singular  (quarter-point)  elements  at  the  tip  of  the  crack. 

However,  for  the  scattering  problem  one  is  generally  interested  in  the  scattered  field  away 
from  this  region.  For  this  purpose  is  it  still  important  to  capture  the  crack-tip 
singularity?  This  paper  addresses  this  question. 
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FORMULATION  AND  SOLUTION 


Consider  a  homogeneous,  isotropic,  and  linearly  elastic  medium  with  a  surface  breaking 
crack  of  arbitrary  orientation  and  shape  as  shown  in  Fig.  I.  This  figure  shows  the 
cross-section  of  the  crack  in  the  xz-plane.  The  stress-free  crack  occupies  the  region 

(x.z)  e  C’.  -  oo  <  y  <  oo.  The  displacement  field  u  (x,y.z;t)  at  any  point  (x.y.z)  and  time 

t  is  taken  to  be  u  (x.zle'fy"'*^*.  where  w  is  the  circular  frequency  and  2ir/f  is  the 

wavelength  in  the  y-direction.  In  this  paper  incident  plane  body  waves  will  be  assumed 
to  be  propagating  in  a  vertical  plane  that  makes  an  angle  0  with  the  x-axis.  Thus  ^  will 
be  determined  by  the  type  of  wave  (longitudinal  or  shear)  and  the  angle  its  direction  of 
propagation  make  with  the  negative  z-axis.  This  angle  will  be  noted  by  6. 

The  problem  of  scattering  of  incident  body  waves  by  the  crack  will  be  solved  by  a 
hybrid  method  which  combines  the  advantages  of  the  finite  element  technique  and  the 
boundary  integral  method.  For  this  purpose  we  consider  two  artificial  boundaries  C  and 
B  (Fig,  2).  The  medium  is  now  divided  into  two  regions.  The  interior  region  R,  is 

bounded  by  B,  part  of  the  free  surface,  and  C*.  The  exterior  region  R,,  is  bounded  by 
the  free  surface  and  C  and  extends  to  infinity  in  the  x  and  z  directions.  The  area 
between  C  and  B  is  shared  by  both  regions. 


Fig.  I.  Geometry  of  the  surface  breaking  crack  showing 
also  the  direction  of  the  incident  field. 


In  regions  R<,  and  R,  the  governing  equation  of  elastic  motion  is  written  as 

Tjj  j  pw*Uj  -  0  (i.j  -  1.2.3)  (1) 

where  Tjj  is  the  stress  tensor,  p  the  mass  density  and  U|  the  ith  displacement  component. 
Solution  to  (I)  satisfying  the  stress-free  boundary  conditions  along  the  surface  of  the  half- 
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space  and  the  crack  surface  is  sought. 


z 


Fig.  2.  Finite  element  grids  around  the  crack. 


Solution  in  the  Exterior  Region  R<i 

In  this  region  the  dispacement  is  composed  of  two  parts. 

Uj  -  uj°^  +  u|*^  (2) 

where  u,(°'  (i  -  1.2.3)  represents  the  free  field  displacement  components  (the  incident 
field  and  its  reflection  from  the  free  surface)  and  u;^**  is  the  scattered  field.  The 
scattered  displacement  field  in  R<,  is  represented  by  a  surface  integral  (Khair  et  aL.  1988), 
after  dropping  the  factor  e'fy. 


Ui(s)(x\  z’) 


C 


(Gij  Tjk-uj  Zjj^)nkdc 


(3) 


Here  Gjj  (x.z;  x'.z'le’fy  is  the  Green's  displacement  tensor  for  the  half-space.  Z|.|^  e'fy  is 
the  corresponding  stress  tensor,  and  n^  defines  the  components  of  the  outv^ard  unit 
normal  vector  to  C.  The  integration  along  C  is  carried  out  in  the  clockwise  direction. 
Expressions  for  G  and  Z  S  have  been  derived  before  (Khair  et  al.,  1988). 
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The  Interior  Region  R, 

This  region  is  divided  into  finite  elements  having  Nj  number  of  interior  nodes  and  N3 
number  of  boundary  nodes.  For  the  finite  element  representation  in  region  R[  the  energy 
functional  is  taken  to  be 
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where  *  denotes  complex  conjugate  and  T.  «  are  stress  and  strain  vectors  defined  as. 
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where  superscript  T  denotes  transpose,  t  and  u  denote  the  traction  and  displacement 

on  contour  B.  respectively.  It  is  assumed  that  the  displacement  field  within  an  element  is 
represented  in  terms  of  the  shape  functions  d|  (x.z)  and  elemental  nodel  displacements 

u^®*  as 
i;l 


11 

Uj^®)  -  y  <fif  u^*J  (i  -  1.2. 
1-1 
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(7) 


The  number  of  nodes  in  each  element  is  given  by  n.  Tjj^®^  and  fjj^®)  are  computed  by 
substituting  (7)  into  strain-displacement  relations  and  these,  in  turn,  into  the  stress-strain 
relations.  Substituting  these  in  (4)  and  taking  variation,  the  equation  of  motion  for 
region  R|  can  be  written  as 
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Elemental  impedance  matrix  [5]^  is  given  by 


(8) 
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Y  is  the  nodal  force  vector  due  to  surface  tractions  on  the  boundary. 

'^B 

Using  the  top  set  of  equations  in  (8)  we  get 
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Matrices  [Bg]  and  [D]  have  been  derived  before  (Khair  et  al.  1988).  Note  that  in  this 
derivation  use  has  been  made  of  the  equation  u  (x.y.z)  -  u(x.z)e*fy. 


Combining  (3)  and  (11)  the  boundary  displacement  U  is  found  to  be 
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Once  tbi  )  IS  found  by  solving  (11).  iU  )  is  found  from  (10). 
B  ~l 
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In  ihe  following  we  present  the  results  for  the  COD  and  z-component  of  the  surface 
displacement  for  incident  longitudinal  (P)  and  vertically  polarized  shear  (SV)  waves.  In 
all  the  numerical  results  presented  it  is  assumed  that  f  •  0.  i.e..  the  waves  are  traveling 
in  the  xz-plaoe.  Also,  to  keep  Ihe  problem  simple  the  crack  is  assumed  to  be  planar  and 
normal  to  the  free  surface.  Results  for  cracks  of  different  orientations  and  for  non- 
planar  cracks  will  be  presented  elsewhere. 


NUMERICAL  RESULTS  AND  DISCUSSION 

The  method  described  above  was  used  to  solve  body  wave  (P  and  SV)  scattering  by  a 
normal  surface-breaking  crack.  The  attention  has  been  focused  on  the  near-field 
displacements  as  they  depend  on  near-the-crack-tip  singular  stress  distribution.  To 
examine  this  dependence  the  region  around  the  crack  was  divided  into  finite  elements 
With  or  without  using  the  quarter-point  singular  elements  near  the  crack  tip.  The 
Poisson’s  ratio  of  the  medium  is  taken  to  be  y  «  1/3.  Displacements  at  all  the  interior 
nodes  were  calculated  for  different  angles  of  incident  using  (II).  Crack  opening 
displacements 


Aj-|u*-u'|  (12) 

Ax-l<-u;|  (13) 

were  calculated.  These  are  shown  in  Figs.  3-6  for  different  frequencies  and  for  angle  of 
incidence  45°.  The  lines  denote  the  results  obtained  by  using  the  crack-tip  singular 
elements  and  t.he  points  without  the  use  of  these  elements.  As  seen  from  these  figures 
the  COD’S  do  not  depend  on  the  accurate  representation  of  the  crack-tip  singular  stress 
field.  That  the  crack-tip  field  does  not  significantly  influence  the  scattered  field  is  also 
seen  from  the  surface  displacement  field  presented  in  Figs.  7  and  8.  Here  the  amplitudes 
of  the  vertical  surface  displacements  near  the  mouth  of  the  crack  are  presented. 

Figures  3-6  also  show  that  at  long  wavelengths  the  COD's  are  nearly  constant  along  most 
of  the  crack  length  changing  rapidly  only  near  the  lip.  This  change  is  influenced  by  the 
way  the  crack-lip  field  is  represented  in  the  finite  element  shape  functions  in  this  region. 
However,  as  noted  above,  this  local  detailed  field  does  not  influence  the  scattered  field 
much,  even  near  the  mouth  of  the  crack. 
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CRACK  OPENING  DISPLACEMENT 


ig-  3,  Distribution  of  along  the  crack  length  for  P  wave 
incident  at  45“  with  the  vertical.  Lines  show  results 

using  crack-tip  quarter-point  signular  elements  and  the 
points  without  that. 


CRACK  OPENING  DISPLACEMENT 
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Fig.  4.  along  the  crack-length  for  the  same  case  as  in  Fig,  3. 


Fig.  6.  ^2  along  the  crack-length  Tor  the  same  case  as  in  Fig.  S. 
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SCATTERED  SURFACE  DISPLACEMENT 
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Normalized  scattered  vertical  displacement  amplitude  near  the 
mouth  of  the  crack  for  P  wave  incident  at  45®  with  ihe 
vertical.  Lines  show  results  obtained  using  the  crack-tip 
elements  and  points  without  them. 


SCATTERED  SURFACE  DISPLACEMENT 
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Same  as  in  Fig.  7  for  SV  wave  incident  at  45"  with  the 
vertical. 
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